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ABSTRACT 



We propose that the unstable periodic variation of the extreme hehum hot 
subdwarf star KIC 10449976 is caused by the stochastic evaporation process 
of a planet orbiting the star with a period of 3.9 days. KIC 10449976 shows 
evidence for photometric periodic modulation of P = 3.9 day. The period is not 
stable, neither in time nor in amplitude. In our model the UV radiation from 
KIC 10449976, whose effective temperature is T^s ~ 40000K, heats the planet 
and causes its envelope to swell and lose mass where dust might be formed. The 
estimated mass outflow rate is ~ 10~^Mjupiter yr~^. Stochastic variations in the 
outflow rate and in dust formation rate lead to the variation in the amount of 
reflected light, hence to the unstable periodic behavior. Self-shielding of the dust 
and the behavior on the day-night boundary might lead to the stochastic outflow 
rate. We predict the presence of a sub-stellar companion to KIC 10449976 with 
a mass of Mjupiter ^ fnp < 20Mjupiter at an orbital separation of ap = 8.3 AU. 
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1. Introduction 



O and B subdwarfs (sdO, sdB) are hot and compact stars that burn helium in their core 
and have a thin envelop e. They are identified with Extreme Horizontal Branch (EHB) stars 



(e.g. 



Heberlll986l 120091 ). with stars evolving off-the zero age HB (Z AHB) through a 'sh ort- 



cut' from the EHB to the WD cooling track (so c alled AGB-Manaue: iDorrnan et al.lll993r). or 



with the merger of tw o white dwarfs (WDs; e.g.. IWebbinklll984j : iHan et al. 



2003; 



Hebei 



2008 



2009; 



Nelemansll2010[ l . The later scenario is attractive when the star is hydr ogen poor. Most 



popular models for the formation of sdBO involve binary interaction (e.g., iHan et al. 



2003 
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Charpinet et al.ll2008l : iGeier et al.ll2010l). These models are supported by the l arge fraction 
of sdB stars in binary systems (e.g.. iMaxted et al.ll200ll : iNapiwotzki et al.ll2004j ). In the first 
two scenarios, those of a star on the HB or evolving off the HE, t he companion is required to 
remove mass from the red giant branch (RGB) progenitor (e.g.. lHan et al.ll2007l ). Many of 
the binary systems have close orbits, implying that the system has gone t hrough a commo n 
envelope (CE) phase during the RGB phase of the sdBO progenitor (e.g., iHan et al.l 120021) . 
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In the case of a CE evolution of a sub-stellar companion to the RGB sta. r , the question 
is whether the sub-stellar compa n ion survived the CE phase (IBear et al.ll201ll : iBear fc Soker 
2OI2I : iTutukoy fc Fedorovall2013r). iBear fc Sokeii (120121 ) suggest that the two Earth-size plan- 



ets found by ICharpinet et al.l ( 1201 ll ) around the sdB star KIC 05807616 are remnant of the 
tidally destructed metallic core of a massive planet. Massive planets might even survive 
the CE phase, as in the case of the planet orbiting the me tal-poor red HB star HIP 13044 
(CD-36 1052) with an orbital period of P = 16.2 ±0.3 days (ISetiawan et al.ll2010r). sdB stars 



with substellar companions at separations of af > 1 AU have also been found ( ISilvotti et al. 



2007; 


Lee et al. 


2009; 


Qian et al. 


2009, 


2012) 



20091 12OI2I ) . In these systems the wide substellar compan- 
ions might hint on a closer planet that went through the CE phase and ejected the envelope 
of the RGB stellar progenitor of the sdOB star. This c loser planet might have been com- 
pletely destructed in the CE process. ISchuh et al.l (120101 ) noted that the increasing number 
of substellar companions to sdB stars may indicate the existence of an undiscovered planet 
population. 

In this paper we examine whether the photometric periodic variation of the sdBO star 
KIC 10449976 can be accounted for by a planet companion. It is an extreme helium- 
rich subdwarf that shows evide nce for photometric periodic modulation of P = 3.9 day 
with an amplitude of ~ 0.02% djefferv et all 120131 ). KIC 10449976 e ffective temp e rature 
is Teff = 40000 ± 300°K and its surface gravity is log^ = 5.3 ± 0.1. I.Tefferv et al.l boi?h 
argue that the sinusoidal modulation is probably an astrophysical variation and not of an 
instrumental origin, and consider possible explanations as follows. (1) Pulsation. The prob- 
lem with pulsation is that for a star with the dimensions of KIC 10449976 the period of 
pulsations sh o uld b e of the order of 200 s <C 3.9 day. (2) Reflection from a close companion. 
Jeffery et al.l ( 120131 ) could find a companion parameters space that can explain the photo- 
metric amplitude and period, as well as the constrain that the radial velocity variability be 
< 50 ± 20 km s~^. However, the lack of stability in the apparent period rule out a stable 
companion. (3) Stellar spots. The photometric variation by stellar spots cannot be rule-out. 
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but the lack of stability in the apparent period is problematic. 



In section[2]we examine the possibility that reflection from an evaporating close-in planet 
(or a low-mass brown dwarf; BD) is behind the photometric modulation of KIC 10449976. 
The planet is heated up, inflates, loses mass, and the reflection from the dust causes the 
variation observed in th e light curve. Evaporation is well studied for pl anets orbiting 



main sequence stars (e.g.. lRappaport et al 



2012 



Perez-Becker fc Chiangll2013l and references 



within), and we use this process to our purposes. The two best known o bserved planets with 
evaporation processes that inf late their atmosphere are HP 20 9458b fIVidal-Madjar et al. 
20031 . l2004h and HD 188733b Jlecavelier Pes Etangs et ahlboioh . The estimated mass-loss 



rates in these cases is 10 g s (lOwen &: Jacksod l2012l ). In section [3] we discuss the 



possible origin of the planet. Our short summary is in section HI 



REFLECTION FROM AN EVAPORATING PLANET 



The heating of the planet by the strong UV radiation from the sdBO parent star can 
cause the planet to inflate and lose mass. The metal rich material in the inflated envelope 
and wind is likely to rapidly form dust. We examine the evaporation with the idea that the 
dust reflects the light from the parent star and account for the photometric variation of an 
un-stable orbital period of 3.9 day. 

We use the term planet here for gas giant planets and low mass brown dwarfs (BD) 
in the total mass range of IMjupiter ^ Mp < 30Mjupiter- The lower limit comes from the 
requirement that the planet survive the common envelope (CE) phase. However, we note 
that lower planet masses are possible if the planet survived the CE phase but already has lost 
most of its mass in the evaporation process. The upper limit comes from the requirement 
that the evaporation rate be high enough to cause the variation in the periodic behavior 
(see below); higher mass BD will have lower evaporation rate. This upper limit is highly 
uncertain, up to a factor of 2. 



T he evaporation of planets has b e en studied before fe.g..lOwen fc Jacksonll2012l : iBear fc Soker 

2011bl : IPerez-Becker fc Chiang|l2013[ l. lOwen fc JacksonI k^vE consider planets around main 
sequence stars and show that planets with orbital separation of a < 0.1 AU will have a hydro- 
dynamical mass-loss due to X-ray and extreme UV (EUV) radiation. The two mechanisms 
that can cause the evap oration are Jean' s escap e (ballistic process) or hydrodynamical pro- 
cess. From the work of lOwen fc JacksonI ( 20121 ) it is evident that the major mechanism for 
mass lose in close planets is hydrodynami cal wind. We below will use an expression for 
evaporation based on our previous study (IBear fc Sokerl l2011bl ). and note that using the 
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expression from lOwen fc JacksonI (|2012| ) gives similar results. 



The effective temperature and sur face gravity of KIC 10449976 are T^s = 40000 ±300X 
and log (7 = 5.3 ± 0.1, respectively ( jJeffery et al.l l2013[ l. These give for the radius and 
luminosity of KIC 10449976 



^Kic = 0.26i?o 



Ml 



KIC 



L 



KIC 



156 



M, 



KIC 



Lcr 



05 



respectively, where Mkic is the mass of KIC 10449976. The orbital separation of the assumed 
planet is calculated from the P = 3.9 day period of the photometric variation 



8.3 



M 



O.5M0 



P 



3.9 day 



P 



0- 



(2) 



Assuming a total light absorption by the planet and emission as a black body from its 
entire surface, the surface temperature of the planet is 



^KIC 



1/2 



KIC" 



5000 



M, 



KIC 



KIC 



(2ap)i/2 \0.5MqJ V40000°K 

where we used the period P = 3.9 day for the planet orbital period 



K 



(3) 



The planet's surface temperature and the strong UV ra diation from the parent star imply 
that evaporation occurs. We use our ( iBear fc Sokeiil2011bl ) previously derived expression for 
evaporation of planets around sdOB stars (eq. 6 and Fig 2 there) and find the evaporation 
rate from the assumed planet around KIC 1044997 6 to be Mp ~ 1 0^^ g s ~^ (the mass-loss 



rate is defined positively). The mass loss rate in iBear fc Sokerl fl2011bl ) is derived for a 



hydrogen rich planet of a mass of Mp = 15Mj, and a radius of Rp = O.IRq, that orbits 
an EHB star having radius, mass, and effective temperature of -Rehb = O.IARq, Mehb = 
O.5i?0, and Tehb = 35500 K, respectively. The ev aporated gas is ionized and can become 
a source of Balmer lines. lOwen fc JacksonI (120121 ) get a similar range of mass-loss rates 
Mp ~ 5 X 10^^ — 10^^ g s~^ for an orbit al separation of 0. 1 AU, and a planet mass of 



Mp ~ IMj (for more details see Fig 9 in lOwen fc JacksonI 120121 ). 



The optical depth of the evaporative flow is given approximately as r 
is the opacity, and 



Krpe, where k 



Pe = 1.8 X 10 



-13 



1014 g 



0.3Rq 



-2 



1 km s^i 



g cm 



-3 



(4) 



is the density at radius r from the planet center. For the evaporative velocity we take 
V = 1 km s~^ which is similar to the value used by iBear fc Soked (120121 ). For the typical 
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radius where most reflection occurs we scale the radius to be up to three times the planet 
radius, r < 3Rp ~ O.Si?©; the radius changes stochastically in according to our model. For 
radiation from the central star we take the typical wavelength to be A ~ lOOOA, and use the 
results of iDraine fc Led (119841 ) for the ISM dust to find the opacity to be /t ^ 10^ cm^ g~^. 
The optical depth if dust forms is then 



1014 g s-1 / V O.3i?0 



1 km s" 



(5) 



Owen &: JacksonI fl2012l ) assume v = 20 km s~i, for r ~ lO^" cm and therefore their 
optical depth is of the order of r ~ 0.4 (Fig. 7 there). For such lower optical depths of 
r > 0.1 the contribution from a more extended region is important, and over all the reflection 
by the evaporative flow is stil signiflcant. Our model can work even for evaporative flows 
speeds of up to ~ 20 km s~i. In any case, th e wind blocks the stella r radiation, and we 
expect the outflow speed to be lower than what lOwen fc JacksonI (120121 ) take. 



The high optical depth, if the typical parameters used here hold, implies several things. 
First, if dust forms it reflects some of the stellar light on an area larger than that of the 
planet. Second, the dust might shield the surface of the star, hence reducing outflow. This 
in turn will block less light, and the temperature of the surface will rise again. A stochastic 
process of outflow might take place. Third, the high temperature calculated in equation (jH]), 
is much above dust formation temperature. This high temperature implies that in order to 
form dust, there must be shielding or rapid expansion of the outflow which cools the gas to 
form dust. Another possibili t y is t hat dust is dragged from the planet by the outflow as 
suggested by lRappaport et al.l (120121 ) to occur in a Mercury-like planet. It is not clear that 
a Jupiter-like planet has enough dust in its outer layers. Any of these processes further 
adds to the stochastic behavior. In the shielding scenario, either dust starts to form and 
shields itself for further dust formation, or dust is formed in the night side of the planet. In 
the later as the planet orbit the star the dust will drag behind and reflects light. Again, the 
reflection percentage is expected to be highly time-dependent. 

Based on the values assumed and derived above, of the mass outflow rate and a size of 
the optically thick dust region of ~ few x O.IRq, we can summarize the main points of this 
section as follows. When partial reflection and a not too extreme inclination of the orbit are 
considered, the dusty outflow derived above can account for an amplitude of up to ~ 0.02%. 
The stochastic behavior of the outflow and dust formation accounts to the n on-stable periodic 



behavior and to the variation of the photometric amplitude. We note that iRappaport et al. 



( 120121 ) suggested that the variable eclipse depths of the K dwarf KIC 12557548 is caused by 
obscuration of dusty evaporative flow from a mercury-like planet. In their model the dusty 
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outflow is variable, hence its optical depth. In our model the variability is in reflection, 
not obscuration, but there are some similarities in the scenario. The largest uncertainty 
in our model is the dust formation process. To better model this process, a numerical 3D 
time dependent numerical code is required that include radiative transfer, chemistry, and 
hydrodynamics. 



3. ORIGIN OF THE PLANET 



There are two scenarios for the evolution of KIC 10449976. The flrst scenario is that 
of a star that lost most of its mass on the RGB and became a blue HB star. Such a star 



evolv es from the HB directly toward the WD cooling curve (AGB Manque; iDorman et al. 



19931 ). and it crosses the region bordering sdB and sdO stars. To be located at the present 
place of KIC 10449976 on t he HR-diagrarn, a hydrogen rich envelope should be of a mass 



O.OO5-O.OIM0 (e.g.,|DormanelaL 



1993 



0stensen 



2009 



2OIOI ) . There are two main 



questions related to this scenario: (1) How did the planet survive the CE phase during the 
end of the RGB phase of the progenitor? (2) Why is the envelope hydrogen poor? 

If our proposed planet is real, it is not the flrs t plane t to raise the question of survivabil- 
ity during a CE envelope phase. ISetiawan et al.l (120101 ) announced the flnding of a planet 
orbiting the metal-poor red HB star HIP 13044 with an orbital period of P = 16.2 ±0.3 days. 
Such a red HB star has an envelope mass of 0.3M(7), which mak es the survivability of the 
planet during he CE phase more puzzle even. iBear et al.l ( 1201 ll ) proposed that the planet 
orbiting HIP 13044 was engulfed by the star as a result of the core helium flash that caused 
this metal-poor star to swell by a factor of ~ 3 — 4. The evolution following the core he- 
lium flash is very rapid, and some of the envelope is lost due to the interaction with the 
planet, and the rest of the envelope shrinks within about a hundred years. This is about 
equal to the spiraling-in time, and the planet sur vived. Mass-loss can increas e the survivabil- 
ity of the planet, as discussed for HIP 13044 bv iTutukov fc Fedoroval (120131 ). In the case of 
KIC 10449976 the envelope mass is much lower, and more likely is that by the time the planet 
spiraled-in to its present orbital separation of ~ 8.3 AU it has expelled most of the envelope; 
the rest of the envelope gas collapsed on the core. Hence, in the case of KIC 10449976 the 
survivability of the planet is easier to explain than in the case of HIP 13044. What might 
further help is another inner planet that collided with the core and helped in expelling the 



RGB envelope (IBear fc Soker 



2011a 



Bear et al.ll201lh 



The helium-rich atmosphere of KIC 10449976 can be related to an inner planet that 
collided with the core, or to the angular momentum of the left-over envelope that collapses 
on the core. Both processes lead to a shear on the core that can dredge-up helium gas to 
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the surface and mix hydrogen to the helium layer . The same processes, including a planet 
companion, were discussed by iDe Marco fc Sokerl (120021 ) to explain hydrogen poor (Wolf- 
Rayet) central stars of planetary nebulae. As well, the mixing o f the hydrogen to the core 
can cause more efficient burning of hydrogen (jSiess fc Lividll999l ). hence reducing hydrogen 
abundance. The influence of rotation on the post-HB evolution requires a detailed study. 

The second scena rio for the for r aation of KIC 10449976 invol yes the merger of t wo H e 



WDs, as discussed by Ijefferv et al.l ( 120131 ) based on the results of IZhang &: Jeffervl ( 



Perets 



2012h. 



2010h. 



In that case t he planet might be formed as a second generation planet (e.g. 
Silvottil (|2008[ ) discuss, as one of the possible less likely evolutionary routes, the formation of 
the sdB star V391 Peg via the merger of two He WDs, and the formation of its planet from 
a disk formed by the r nerger process. In the c ase of the two Earth-sized bodies orbiting the 
sdB star KIC 05807616 lCharpinet et al.l (120111 ) considered the formation of second generation 
planets after the merger of two WDs to be an unlikely channel. 



4. SUMMARY 



In a recent paper (jjeffery et al.ll2013[ ) report a 3.9 day periodic photometric variation in 



the hot subd warf (sdBO) s t ar KI C 10449976. The period is not stable, neither in time nor 



in amplitude. iJefferv et al.l ( l2013l ) discuss three possible scenarios as possible explanations: 
stelar spots, pulsation, and their favorable explanation of reflection from a stellar companion. 
However, they don't account for the unstable periodic behavior. 

In this paper we offer an alternative scenario for the unstable periodic behavior of 
KIC 10449976 based on a planet (or a low-mass brown dwarf) that orbits KIC 10449976 at 
an orbital separation of = 8.3 AU. We argue that the UV radiation from KIC 10449976, 
whose effective temperature is Tgg = 40000 ± 300K, heats the planet and causes its envelope 
to swell and lose mass. We assume that dust is formed in the outflowing gas. This gas reflects 
the light from KIC 10449976 and form the photometric periodic variation. Our estimated 
mass-loss rate, section [21 and the assumption of dust formation implies a large enough optical 
depth, equation to account for the required amount of reflection. 

In our model, stochastic variations in dust formation rate lead to the variation in the 
amount of reflected light, and hence to the deviation from a strictly periodic behavior, both 
in period and amplitude. Several effects can lead to the stochastic dust formation process. 
First is the self-shielding of the dust. As more dust is formed it block the stellar surface, 
and evaporation rate, hence dust formation, decreases. This feedback effect can lead to a 
stochastic behavior. A second effect is the day-night variations, in particular if the rotation 
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of the planet is not completely synchronized with its orbit. Dust might be formed due to 
'weather' effects, in particular on the day-night boarder. Dust that is formed from the hot 
region there can leave the planet, as in a comet. 

There are several processes that must be studied in greater detail to put our proposed 
scenario on a firmer ground. First is dust formation from an evaporating planets in such an 
intense UV radiation. The calculation is complicated as the region of day-night boundary 
must be taken into account. Also, the evaporation rate must be better determined. And 
then the amount of reflection should be calculated. 

Our prediction is simply that there is a secondary body there with a mass of Mjupiter ^ 
iT^p ^ 20Mjupitcr at an orbital separation of ap = 8.3 AU. We encourage accurate spec- 
troscopic study of KIC 10449976 in order to look for such a planet. The expected orbital 
velocity of the parent star KIC 10449976 is ~ 0.2 — 4 km s^^. As well, some weak Ha 
emission from the outflowing gas might be detected, if the planet's atmosphere is hydrogen 
rich. The Ha emission will be very weak, much lower than the reflected light. 

This research was supported by the Asher Fund for Space Research at the Technion, 
The US - Israel Binational Science Foundation. 
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